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Abstract 

(D ■ 

C/^ I We provide an equation of state for high density supernova matter by applying a momentum- 

dependent effective interaction. We focus on the study of the equation of state of high-density 
. and high-temperature nuclear matter containing leptons (electrons and neutrinos) under the 

chemical equilibrium condition. Thermal effects on the properties and equation of state of 
nuclear matter are evaluated and analyzed in the framework of the proposed effective inter- 
action model. Since supernova matter is characterized by a constant entropy we also present 
the thermodynamic properties for the isentropic case [1]. 



1 Introduction 



> 

^ ■ Knowledge of the properties of the equation of state (EOS) of hot asymmetric nuclear matter is 
! of fundamental importance to understand the physical mechanism of the iron core collapse of a 
^ massive star which produces a type-II supernova, and the rapid cooling of a new born hot neutron 
0^ star. Additionally, the EOS defines the chemical composition, both qualitative and quantitative, 
^ of the hot nuclear matter [21 El HI [5] . Supernova matter which exists in a collapsing supernova core 
O ■ and eventually forms a hot neutron star at birth is another form of nuclear matter distinguished in 
>- ! the participation of degenerate neutrinos and electrons [6j. It is characterized by almost constant 
^ ' entropy per baryon S" = 1 — 2 (in units of the Boltzmann constant ks) throughout the density n 
^ [ and also by a high and almost constant lepton fraction Yj = 0.3 — 0.4 in contrast with ordinary 
neutron star matter where S = and Yj < 0.05. These characteristics are caused by the effects of 
neutrino-trapping which occurs in the dense supernova core where a neutron star is formed. 

This work is a continuation of our previous work concerning the EOS of hot /3-stable nuclear 
matter in cases where neutrinos have left the system [7j. More specifically, in order to study 
the properties and the EOS of hot nuclear matter, a momentum-dependent effective interaction 
model (MDIM) has been applied, one which is able to reproduce the results of more microscopic 
calculations of dense matter at zero temperature and which can be extended to finite temperature 
[31 [71 El [9]. The main incentive for the present study is the fact that only few calculations of 
the equation of state of the supernova matter at high densities are available, although at lower 
densities (n < no) (where no is the saturation density) reliable results are already available. 

2 The model 

The model we use here, which has already been presented and analyzed in our previous papers 
[3, [HI [HI [lO], is designed to reproduce the results of the microscopic calculations of both nuclear 



1 



and neutron-rich matter at zero temperature and can be extended to finite temperature [31 II]. 
The energy density of the asymmetric nuclear matter (ANM) is given by the relation 



e(n„, Up, T) = efei„(n„, T) + el^^iup, T) + Vint{nn, Up, T), 
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where (n^) is the neutron (proton) density and the total baryon density is = n„ + Up. The 
contribution of the kinetic parts are 



(27r)3 2m 



(/n(n„, k, T) + fp{np, k, T)) 



(2) 



where fr, (for t = n,p) is the Fermi-Dirac distribution function. 

Including the effect of finite-range forces between nucleons, the potential contribution is pa- 
rameterized as follows [3] 
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g{k,Ai)frinr,k,T). 
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In Eq. ([3]), / is the asymmetry parameter (/ = (r2„ — np)/n) and u = n/nQ, with Hq denoting 
the equilibrium symmetric nuclear matter density, uq = 0.16 fm~'^. The asymmetry parameter / 
is related to the proton fraction Yp by the equation / = (1 — 2Yp). The parameters A, B, a, Ci, C2 
and B' which appear in the description of symmetric nuclear matter are determined in order that 
E{n = tiq) — mc^ = —16 MeV, Uq = 0.16 fm~^, and the incompressibility are K = 240 MeV. The 
additional parameters Xq, X3, Zi, and Z2 used to determine the properties of asymmetric nuclear 
matter are treated as parameters constrained by empirical knowledge 

The function, g{k,Ai), suitably chosen to simulate finite range effects, has the following form 



g{k,A, 



' k_ 
A" 



(5) 



where the finite range parameters are Ai = l.bk^p and A2 = 3^^^ and k^p is the Fermi momentum 
at the saturation point tt-q. 

The energy density of asymmetric nuclear matter at density n and temperature T, in a good 
approximation, is expressed as 



where 



e(n, T, I) = e{n, T, / = 0) + esymin, T, I), 



, T, /) = nl'E%l{n, T) = nl' {E^;i{n, T) + Ef^{n, T)) . 
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In Eq. (I7j) the nuclear symmetry energy E^°y^{n,T) is separated into two parts corresponding to 



the kinetic contribution E^y^{n,T) and the interaction contribution E^^l^{n,T). 
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From Eqs. ^ and ([7]) and setting / = 1, we find that the nuclear symmetry energy E^°y.^{n, T) 



is given by 
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[n, T) = - (e(n, T, / = 1) - e(n, T, / = 0)) . 
n 



Thus, from Eq. and by a suitable choice of the parameters Xq, x^, Zi and Z2, we can obtain 
different forms for the density dependence of the symmetry energy El°l^{n,T). 

In a very recent work, QjJ the authors carried out a systematic analysis of the nuclear symmetry 
energy in the formalism of the relativistic Dirac-Brueckner-Hartree-Fock approach. In this case 
Esymiu) is obtained with the simple parametrization Egymiu) = Cu' with 7 = 0.7 — 1.0 and 
C ~ 32 MeV. The authors concluded that a value of 7 close to 0.8 gives a reasonable description 
of their predictions although the use of different functions in different density regions may be best 
for an optimal fit pJJ. The results of Ref. [llj are well reproduced by parameterizing the nuclear 
symmetry energy according to the following formula 

E*°*^(n,T = 0) = 13«2/3 + 17F(«), (9) 

where the first term of the right part of Eq. ^ corresponds to the contribution of the kinetic 
energy and the second term to the contribution of the interaction energy. 

For the function F{u), which parameterizes the interaction part of the symmetry energy, we 
apply the following form 

F{u)=u. (10) 

The parameters Xq, X3, Zi and Z2 are chosen so that Eq. ([H]), for T = reproduces the results of 
Eq. ([9]) for the function F{u) = u. 



2.1 Thermodynamic description of hot nuclear matter 

In order to study the properties of nuclear matter at finite temperature, we need to introduce the 
Helmholtz free energy F which is written as 



F(n, T, /) = E(n, T, J) - TS{n, T, /). (11) 

In Eq. (11 II) . E is the internal energy per particle, E = e/n, and S is the entropy per particle, 
S = s/n. From Eq. flTTl) it is also concluded that for T = 0, the free energy F and the internal 
energy E coincide. 

The entropy density s has the same functional form as that of a non-interacting gas system, 
given by the equation 

/d h 
^ [fr In /. + (1 - /.) ln(l - fr)] , (12) 

while the pressure and the chemical potentials defined as follows [12] 

\ / S,Ni \"'*V5y,n,^, 

At this point we shall examine the properties and the EOS of nuclear matter by considering 
an isothermal process. In this case, the pressure and the chemical potentials are related to the 
derivative of the total free energy density / = Fn. More specifically, they are defined as follows 

It is easy to demonstrate by applying Eq. ( fT4l) that (see for a proof [I3] as well as [14] ) 



P/ n,T 
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A = ^"-^- = -(^)^^^- (15) 

We can define the symmetry free energy per particle Fsym{n,T) by the following parabolic 
approximation (see also p3l [15] ) 

F(n,T,/) = F(n,T,/ = 0) + /2F,j,„(n,T) 

= F(n,r,J = 0) + (l-2Fp)2F,,^(n,T), (16) 

where 

i^syn^(n, T) = F(n, T, / = 1) - F(n, T, / = 0). (17) 

It is worth noting that the above approximation is not valid from the beginning, but one needs to 
check the validity of the parabolic law in the present model before using it. In Ref. ^ we have 
proved the validity of the approximation ( fT6l) . 

Now, by applying Eq. ( fT6l) in Eq. ( |T5|) . we obtain the key relation 

/i = /i„ - /ip = 4(1 - 2Yp)F,yUn, T). (18) 

The above equation is similar to that obtained for cold nuclear matter by replacing Esym{n) with 



2.2 equilibrium in hot proto- neutron star and supernova 

Stable high density nuclear matter must be in chemical equilibrium with all type of reactions, 
including the weak interactions in which f3 decay and electron capture take place simultaneously 

n — > p + + Ue, p + — > n + z/g. (19) 

Both types of reactions change the electron per nucleon fraction, Yg and thus affect the equation 
of state. In a previous study, we assumed that neutrinos generated in these reactions left the 
system [7]. The absence of neutrino-trapping has a dramatic effect on the equation of state and 
is the main cause of a significant reduction in the values of the proton fraction Yp |6l [16]. The 
equation of state of hot nuclear matter in /3-equilibrium (considering that it consists of neutrons, 
protons, electrons and neutrinos) can be obtained by calculating the total energy density etot as 
well as the total pressure Ptot- The total energy density is given by 

etot{n, T, I) = et{n, T,I)+ ^ Q(n, T, J), (20) 

where ef,(n, T, J) and ei{n,T,I) are the contributions of baryons and leptons respectively. The 
total pressure is 

PUn,T,I) = P,{n,T,I)+ ^ Pi{n,TJ), (21) 

where Pb{n, T, I) is the contribution of the baryons i.e. 

Pbin, T,I)=T Yl Sr{n, T,I)+ ^ n,/i,(n, T, /) - eb(n, T, I), (22) 

T=p,n T=n,p 

while Pi{n,T,I) is the contribution of the leptons. From Eqs. ([201) and ([2T!) we can construct 
the isothermal curves for energy and pressure and finally derive the isothermal behavior of the 
equation of state of hot nuclear matter under /^-equilibrium. 
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3 Results and Discussions 



We calculate the equation of state of hot asymmetric nuclear matter by applying a momentum 
dependent effective interaction model describing the baryons interaction. We consider that nuclear 
matter contains neutrons, protons, electrons and neutrinos under /3-equilibrium and charge neu- 
trality. The key quantities in our calculations are the proton fraction Yp and also the asymmetry 
free energy defined in Eq. (ITTl) . It is worth pointing out that since the supernova explosion itself 
is a dynamic phenomenon, the chemical composition of matter changes according to the evolution 
of the star all the time |17]. During supernova explosion, the chemical composition of matter 
reaches equilibrium not in the whole star but locally. In our present work we assume matter in 
the chemical equilibrium for simplicity in order to analyze the properties of hot neutron star and 
supernova matter. 

Following the discussion of Takatsuka et al. [6J, we attempt to extend the discussion concerning 
the dependence of equilibrium fraction Ye{= Yp) on the baryon density as well as on the nuclear 
symmetry energy. We ignore the temperature effect to clarify the situation. Actually, the situation 
does not change by including finite temperature effects. The energy per baryon of supernova matter 
EsM and cold neutron star matter Ej^s are expressed as function of n and Yp (see also ref. [6]) as 

EsM{n,Yp) = Eb{n,Yp) + E,{n,Yp) + E,^{n,Yp) (23) 

= E,{n,Yp = 0.5) + EsyUn)il-2Ypy 

+ 253.6M^/Vjf/^ + 319.516m^/3(Y, - 

ENs{n,Yp) = Eb{n,Yp) + E,{n,Yp) (24) 

= E,{n, Yp = 0.5) + EsyMil - 2Yp)^ + 253.6u'/%^/\ 

Esym{f^) is plotted in Fig. 1(a) for the three different parametrizations. In the same figure we have 
included recent results provided in reference [TT] achieved by performing microscopic calculations 
in asymmetric nuclear matter. In this case Esym{n) is obtained with the simple parametrization 

with 7 = 0.8 and C = 32 MeV. It is obvious that the results of the above parametrization, 
correspond very well with the parametrization F{u) = u which is proposed here. 

The equilibrium proton fraction Yp is calculated by solving the equation dEsM,Ns/dYp = for 
various values of the density n, Esym{n) and Yj = 0.4 for supernova matter (see Fig. 1(b)). In the 
case of cold neutron star matter, Yp depends strongly on both the baryon density and the values 
of the Egymin). This is not the case for supernova matter where the effect of nuclear symmetry 
energy in determining Yp is less important than in cold neutron star matter. In addition, Yp, for 
a fixed parametrization of F{u) is almost constant with respect to n. 

In Fig. 2 we plot the contribution of the baryons Sb, leptons Si and the total Stot to the entropy 
per baryon. In all cases, S" is a decreasing function of the baryon density n. Temperature affects 
appreciably both baryon and lepton contribution. It should be noted that the contribution of 
baryons Sb may be written as Sb = Skin + Smt, where the term Skin originates from the temperature 
effect on the kinetic part of the energy density and Sint reflects thermal effects on the potential 
energy density. 

In Fig. 3, we display the contribution to internal energy E from baryons Eb and leptons Ei 
for Yi = 0.3 and Yi = 0.4 and for various values of T. The most striking aspect is that the lepton 
energy, Ei = Eg + E^^, dominates in the internal energy of the matter up to n ~ 0.6 fm~^ (for 
Yi = 0.3) and n ~ 0.8 fm"^ (for Yi = 0.4). This is a characteristic of the supernova matter and is 
in remarkable contrast with the situation of cold neutron star matter [B]. The contribution from 
baryon Eb gets larger with the increase of n and is comparable with Ei for high values of n. 
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Figure 1: a) The nuclear symmetry energy for three different parametrization of the interaction 
part with the results of reference [TT] (see text for more details) and b) the proton fraction Yp versus 
baryon density for cold neutron star matter (down curves) and supernova matter (up curves) for 
the three different parametrization of the nuclear symmetry energy. 




Figure 2: Contribution to the total entropy per particle of a) baryons (Sb), b) leptons (Si) and c) 
the total entropy (Stot) versus the baryon density for various values of T for total lepton fraction 
Yi = 0.3 and Yj = 0.4. 




Figure 3: Contribution to the total energy per particle of baryons (Ej,), leptons (£"/) and the total 
energy (Etot) versus the baryon density for various values of T for total lepton fraction a) Yi = 0.3 
and b) y; = 0.4. 
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Figure 4: The pressure of baryons (Pb), leptons (P/) and the total pressure {Pm) versus the baryon 
density for various values of T for the ) Yi = 0.3 and b) Yi = 0.4. 




Figure 5: Temperature T-density n relation with Yi = 0.3 (solid line) and Yi = 0.4 (dashed line) 
for 5" = 1. Stars denote the results for the case with S = 1 and Yi = 0.4, extracted from the 
results by Lattimer et al. ^19j . 
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Figure 6: Internal energy per baryon versus density n for dense supernova matter (SM) in com- 
parison with that of cold neutron star matter (NS) by applying the same model. The case of 
supernova matter corresponds to S* = 1 and Yi = 0.3. The contribution of each species is plotted 
separately. 
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The contributions of baryon and leptons on the total pressure are presented in Fig. 4. In con- 
trast to the situation of the internal energy, the nuclear part contribution plays a more important 
role compared with the lepton part. The lepton pressure Pi is comparable to baryon pressure Pb 
up to 72 ~ 0.2 fm~^, but for higher values of n it is significantly small. 

As pointed out by Bethe et al. [18], the crucial feature in determining the evaluation of a 
collapsing pre-supernova core is that the entropy per particle is very low, of the order of unity (in 
units of the Boltzmann constant fc^), and nearly constant during all the stages of the collapse up 
to the shock wave formation. Therefore, the collapse is an adiabatic process of a highly ordered 
system. So, since the supernova matter is characterized by a constant entropy and constant lepton 
fraction, we shall also discuss the properties under this condition. This can be done by converting 
the results for isothermal case (T=const) into those for adiabatic case (5'=const) in terms of the 
T — n relation constrained by a constant entropy. 

The T = T{n) relation is constructed by {T, n} values to satisfy 5'(n, T)=const in an 5* — n 
diagram. Fig. 5 shows the results for = 0.3 and Yi = 0.4 for S = 1. Temperature is an increasing 
function of n. Furthermore, for the same density, the temperature is higher for lower values of YJ. 
The values of T for various values of n are derived, for the two cases, with the least-squares fit 
method and found to take the general form 

T(n) = an!', 

where a = 35.412, b = 0.70481 for Yi = 0.3 and a = 32.35706, b = 0.67694 for Yi = 0.4. The 
results of this study are in very good agreement with those of Takatsuka et al. [6J. The stars at 
lower density denote the {T, n} values for S = 1 and Yi = 0.4 which are derived from Lattimer et 
al. [19]. It is concluded that the temperature increases considerably when moving from the outer 
part of the star to the center in order to maintain a constant value of the entropy per baryon. 

Finally, in Fig. 6 we compare the EOS's between supernova matter and cold neutron star 
matter. The case for supernova matter corresponds to S' = 1 and Yi = 0.3. It is thus clear that 
the internal energy Etot of supernova matter (SM) is remarkably larger than that of neutron star 
matter (NS). As far as the nucleon part Eb is concerned, the Eb in SM is slightly lower than 
that in NS due to the large energy gain in symmetry energy (see also |6j). However, the lepton 
contribution on the internal energy Ei is remarkably larger in SN matter compared to NS matter 
due to the effect of a large lepton fraction, that is, a large kinetic energy of abundant leptons. 
High temperature also contributes to the stiffening, but it is less effective than the high lepton 
fractions. 



4 Summary 

The evaluation of the equation of state of hot nuclear matter is a major challenge for nuclear physics 
and astrophysics. EOS is the basic ingredient necessary for studying the supernova explosion as 
well as for determining the properties of hot neutron stars. The motive for the present work has 
been to apply a momentum- dependent interaction model for the study of the hot nuclear matter 
EOS under /3-equilibrium. Special attention has been dedicated to the study of the contribution 
of the components of /^-stable nuclear matter on the entropy per particle, a quantity of great 
interest in the study of structure and collapse of supernova. The above EOS can be applied to 
the evaluation of the gross properties of hot neutron stars i.e. mass and radius. 
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